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Abstract 

Vibrating systems with singular mass-inertia matrices arise in recent continuum 
models of Smart Structures (beams with PZT strips) in assessing the damping at- 
tainable with rate feedback. While they do not quite yield “distributed” controls, we 
show that they can provide a fixed nonzero lower bound for the damping coefficient at 
all mode frequencies. The mathematical machinery for modelling the motion involves 
the theory of Semigroups of Operators. We consider a Timoshenko model for torsion 
only — a “smart string,” where the damping coefficient turns out to be a constant at 
all frequencies. We also observe that the damping increases initially with the feedback 
gain but decreases to zero eventually as the gain increases without limit. 
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1. Introduction 


We consider a class of vibrating systems with singular mass-inertia matrices. This 
would appear to be a characteristic feature of continuum models of smart structures 
[1-4] — beams with self-straining material. In this paper in order to keep the exposi- 
tion simple and illustrate the nature of the results, we limit ourselves to a Timoshenko 
model for torsional motion only — a “smart string” — with rate feedback. We show 
that the damping coefficient is a constant at all frequencies. Also the damping in- 
creases as the gain increases initially but decreases to zero eventually as the gain 
increases without limit. The description of the motion dynamics involves the use of 
theory of Semigroups of Operators. 


2. Timoshenko Model: Torsion with Rate Feedback 

We consider a simplified Timoshenko “smart” beam confined to torsion mode only 
— clamped at one end and self- straining material along the entire length with rate 
feedback. The continuum model dynamics are then described by (see [4]): 

p9(t, s) — c9"(t, s) — 0, 

9(t,L) = 0 
C 0'M) + a0(t,O) = 

Here 9(t, s) is the torsion angle and p the inertia parameter and c the stiffness pa- 
rameter; and ct > 0 is the feedback gain. The dots indicate derivative with respect to 
time and the primes the derivative with respect to s, the position along the beam, L 
being the beam length. We may consider this as a “singular” version of the “regular” 
problem where we do have a nonzero “tip-inertia” — with the boundary equations 

m9(t, 0) -I- c9'(t, 0) -I- a9(t, 0) = 0 

considered for example in the original version of the SCOLE [5]. We have 

m = 0 

in the present case. To underscore the “singular” feature, let us consider first how 
the problem is treated for the regular case: 

m ^ 0 . 

As in [5], we begin with the “state” space: 

n = l 2 [o,l] x e 1 


0 < s < L\ 0 < t 


(2.1) 
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including in other words the boundary 


x = 


/(•) 

b 


value in the state. We use the notation: 
/(•) € L 2 [0,L]; b € E l . 


We define next the stiffness operator A by 


V(A) 



/(•) 

X — 

b 


/'(•)./"(•) €L 2 [0,L]; b = / (0) ; f(L) = 0 


and 

- c /"(0 

— 

-c/'(0) 

Thus defined /I is self-adjoint and nonnegative definite and 

r L 

[Ax, x] = cl |/'(s)| 2 ds = Potential Energy. 

We can therefore define the positive square root, y/~A , and we note that 

v(Sa) = 

Define the mass-inertia operator M by 

M x = 


x 


f 

/( 0 ) 


/'(•)eL 2 [0,L]; f(L) = 0 


Pf(-) 

mb 


and note that M is nonsingular, m being nonzero. Define the control operator B on 
E 1 into by 

Bu = 0 

u 

Then the equations (2.1) go over into the abstract version: 

M x(t) + Ax(t) -I- aBB*x(t) = 0. (2-2) 


Note that this corresponds to the “collocated” case, with rate feedback. To proceed 
further we need to construct the energy space 

He = V (\/I) x H 

in which the inner product is the energy inner product: 
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Defining 


y = 


x x 

*2 


z = 


21 

22 


[Y,Z]e = [Va x lt sFA 2ij + [Ml 2. 2 2] 

~ Potential Energy + Kinetic Energy . 


0 

-M~ l A 


I 


(2.3) 


(2.4) 


V(A) = 


' 

fll(-) 

y\ = 

M-) 

€ V(A) 

v 

6,(0) 


I — 

fla(-) 

e 2 (0) 

1 

y 2 = 

02(') 
0 2 ( 0 ) 

€ H(7I) 


w 

02(0) 

AY = 

“ m 0 2(O) 

we see that 

Re[AY, Y) e = -a||B*y 2 || 2 = -a|W)| 2 
and hence A is dissipative and generates a contraction co-semigroup: 


(2.5) 


IIVMIIs < lir(0)|| fi . 

We can show that 

||y(t)IU ^0 as ( -* 00 

but the semigroup is not exponentially damped. The resolvent 7£(A,,4) is compact. 
The eigenalvues of A are the roots of 

D( A) = (a + 77 i A) Sinh XvL + cv Cosh \vh =0; v 1 = ^ . (2.6) 

We note that D( A) is an entire function of order one and completely regular growth 
[6] and hence has a countable number of eigenvalues. Denoting them by {A*} we have 

A k = -|crfc| ± ioJk\ kfcl t 00 
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and [5,6]: 


M 


EM = £ 


— ♦ 0, not monotonic. 


Thus the damping due to feedback goes to zero as mode frequency increases. The 
damping increases as a increases and as m decreases. 


Singular Case 

With this background, let us now get back to our problem where 

m = 0 . 

We see that (2.5) no longer is defined, since m is zero. M is singular and hence we 
can no longer use (2.3) as an inner product. Hence we need to change the definition 
of the energy space He- We note that the total energy is in this case given by 


ds. 


Hence we define now 


Y = 


[ L \e'{t,s)\Us + P [ L \o(t, s )f 

JO Jo 

He = V(/A) x L 2 [0,L] 

e v(/a) 


h(-) 

/i(0) 

/*(•) 


/i(-) 

/i(0) 


y i = 

2/2 = /2(-) € L 2 [0, L] 
[Y,Y] e = Va yi + p[j/2» 1 / 2 ] 


Define A by: 


V(A) = 



/i(-) 

Y = 

/i(0) 


/a(-) 


2/i = 


/i(-) 

/i(0) 


/»(*) 

f/((0) 


€ V{A) 
G v(Va) 


AY = 


f 2 (‘) 
if m 


a 

C fft 


jfl'(-) 


4 



We can then verify that 


Re [AY, Y) = — |/((0)| 2 < 0. 

Q 

Hence A is dissipative, generates a co-semigroup which is a contraction, and (2.1) 
goes over into: 

Y(t) = AY(t) 

and 

limn* < \\y(o)\\e. 

The resolvent TZ(X, A) is again compact. The eigenvalues are the roots of 

D{ A) = cv Cosh XvL + aSinhAj/L = 0 (2.7) 

where we note that D( A) is an entire function of order one and of completely regular 
growth [6] with a countable number of eigenvalues {A*}: 

A it = -|<7*| ± iu k , |w fc | T oo- 


In fact |<Tjt| is given by 


•"‘I = 2 Tv L ° g 


Q T cv 


a — cv 


= M 


and 


Ait = 


— |cr| ^ i(2k + 1) 7T 
Lv 

— \a\ iku 

Lv Lv 


— , a < c v 


Lv 2 

Q > CV. 


The eigenfunctions are 


0 k = A k Sinh X k v(L — s), 0 < s < L 

$it = A k 


U) 
9 k { 0) 

Ait^jt(‘) 


= A;t$it 


*k = A k 


Ok(-) 

Moy 

— Afe ^jfc(-) 


(2.8) 
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A'Vk = a *#*. 


We note that {$ fcl 'f'*} yield a Riesz basis for all values of a, and in turn we have a 
bi-orthogonal modal expansion for any Y in "He- 


y = E 

k 


[**, #fc] 




The semigroup generated by ,4 is exponentially damped: 

l|r(()ll s < l|r(0)|| B e“ w ‘. 

Finally we note the dependence of the damping coefficient on the gain a. We have 
that the damping increases as a increases from zero to cv. At a = cv, we have 

M = 00 

and for a > cv, \a\ decreases as a increases. The function is plotted in Figure 1. The 
apparent paradox of decreasing damping can be explained by noting that as a — * 00 , 
we are in the neighborhood of q = 00 and here the modes are determined by the 
roots of 

SinhAvL = 0 

and for large a we may write: 

cv 

Sinh \vL + — Cosh A vL = 0 

a 

and consider A in place of a. In practice, the attainable gain is such that 

Q 

— < 1 

cv 

These results extend to the more general anisotropic Timoshenko models consid- 
ered in [7] — the damping coefficients are no longer constant but tend to a nonzero 
limit as the mode frequency increases. Details of this analysis will appear elsewhere. 
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Damping vs. Gain 


Damping 

(<tLu) 



GAIN (a/a/) 
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